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ABSTRACT 

The polarization of gamma ray emission scattered off the baryon rich material that collimates a GRB 
fireball and the coasting speed of the irradiated matter are calculated numerically for different geometries 
of the radiation source and the collimating wall. It is shown that when the scattering material is Compton 
sailing, the direction of maximum polarization is quite generally well within the beaming cone of scattered 
radiation. As a result, the probability for observing bright highly polarized GRB's increases considerably, 
provided the Lorentz factor of the coasting matter is not well below 30, and the scattered radiation is 
highly polarized even when the beam thickness is large compared to 1/T. It is suggested that correlation 
between polarization and intensity could provide clues as to whether energy flows from matter to photons 
or the reverse. 

Subject headings: black hole physics — gamma-rays: bursts and theory 



1. INTRODUCTION 

Gamma ray bursts (GRB's) are believed to contain 
about 10 51 ergs in the form of several x 10 57 photons aver- 
aging several hundred KeV in energy. On the other hand, 
they are believed to have Lorentz factors of at least 10 2 , 
and, hence, at most ~ 10 51 protons. For each proton, 
there must be at least 10 5 photons/baryon emitted. The 
question arises: what determines the number of photons 
emitted per baryon and why should the average photon 
energy in a GRB be consistently close to the electron rest 
mass? 

A popular model posits that these photons are emit- 
ted as synchrotron radiation in an optically thin envi- 
ronment, typically at distances 10 13 to 10 15 cm. How- 
ever, this begs the question of why a) the synchrotron 
frequency in the lab frame, which depends on the mag- 
netic field, b) the comoving Lorentz factor 7' of the emit- 
ting electrons, c) the bulk Lorentz factor T, and d) the 
viewing angle together so frequently conspire to put the 
peak luminosity near 300 KeV. The fact that this peak 
energy has been observed to vary from ~ 2 MeV all the 
way down to ~ 10 KeV notwithstanding, we contend that 
the typical value of several hundred KeV is remarkable. 
For example, the strong linear correlation in the HETE II 
data set (Atteia et al. 2003) between the square root of 
the isotropic equivalent fluence Ei SO and the spectral peak 
hv P eak, hu peak /100KeV ~ [E^o/lO^erg] 1 / 2 , suggests that 
the softer GRB's are intrinsically weaker and quite possi- 
bly merely the peripheries of GRB jets. The bright bursts 
show remarkably little variation in the location of the spec- 
tral peak. 

A detailed analysis by Preece et al. (2002) shows that 
even if the spectral peak distribution could be arranged 
within the optically thin synchrotron model, the spectra 
at low energy, which typically show dn^/dE oc E" 1 , are 
inconsistent with optically thin, efficiently cooling syn- 
chrotron emission which would give dn 1 jdE oc E~^l 2 . 
(It could be argued that the bulk Lorentz factor T is fre- 



quently high enough that the synchrotron cooling is inef- 
ficient, which would give a harder spectrum, but then it 
becomes hard to understand how the emission, which de- 
pends on a high power of T in this regime, is frequently so 
efficient.) The recent claim of high polarization in GRB's 
(Coburn and Boggs 2003) has nevertheless strengthened 
the views of some that the emission mechanism is optically 
thin synchrotron (e.g. Lyutikov, Parlev and Blandford 
2003). Others have argued that the claimed degree of po- 
larization, 70 to 80 percent (but with considerable exper- 
imental uncertainty), is too high to be accommodated by 
this mechanism, which limits the polarization to at most 
56 percent for the case of a wide beam, A9 3> 1/T. 

At the time of this writing, the observational situation 
is, in the opinion of the authors, unresolved on several 
fronts. The significance of the X-ray flashes is not yet 
fully understood. The degree of polarization estimated by 
Coburn and Boggs may be considered uncertain by sev- 
eral factors: e.g. large error bars, the chance proximity of 
the modulation frequency of the RHESSI detector and the 
intrinsic variability of GRB021206, including 2s variabil- 
ity in the hardness (Hajdas et al. 2003), and the absence 
of other bursts to date that could establish how typical 
this burst was. Much of the motivation for the model dis- 
cussed in this paper is thus tentative, but it should be 
clear that working out the predictions theoretically is part 
of the process of resolving the issues observationally. 

An alternative model has been developed in a series of 
papers (Levinson and Eichler, 1993, 2000, 2003; Eichler 
and Levinson 1999, 2000, 2003; Eichler 1994, 2003; van 
Putten and Levinson 2002, 2003) originally and still moti- 
vated by the above questions. Its basic assumption is that 
the extreme baryon purity that must exist in the high T 
GRB outflows is enforced by the event horizon of a cen- 
tral black hole, so that the original baryonic content in the 
outflow must vanish. The photosphere near the axis of the 
outflow is therefore at about the annihilation radius of the 
pairs, which until then constitute the main component of 
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the fireball plasma. The fireball's ability to manufacture 
more photons drops precipitously beyond the e ± annihi- 
lation radius, so the total photon entropy of the GRB is 
frozen near its value at pair annihilation, when the average 
energy per photon is of the order of magnitude of m e c 2 . 
[Actually, the comoving temperature should be about an 
order of magnitude (the log of the compactness parame- 
ter) less than m e c 2 if there is perfect thermal equilibrium, 
but the dissipation introduced by the friction of the fire- 
ball with the surrounding walls probably introduces some 
non-equilibrium signature on the spectrum, and the bulk 
Lorentz factor is probably at least a few even in the pres- 
ence of such friction.] 

In this class of models, the average photon energy in a 
GRB, though somewhat dependent upon the geometry of 
the walls and how the photons bounce off them, is natu- 
rally of order m e c 2 without fine tuning of the parameters. 
As long as photons can bounce off the walls at large angles 
and make pairs with other outgoing photons, the dissipa- 
tion continues, the expansion is accompanied by additional 
photon production, and the bulk Lorentz factor remains 
moderate. As soon as the energy density is below that 
needed to maintain pairs, transparency sets in and the 
photon entropy remains nearly constant thereafter. 

Yet another class of models posits that GRB's are pho- 
tons that are upscattered optical/UV photons by relativis- 
tic electrons with a Lorentz factor of order 300 (Shaviv 
and Dar 1995a, b; Lazzati et al. 2000; Dado, Dar, and De 
Rujula 2003). In these models, the number of photons is 
conserved and the energy of the outgoing photons, which is 
provided by the Compton scattering material, has been in- 
creased by a factor of 10 5 to 10 6 . Thus, the once scattered 
photons constitute the primary GRB beam. Because the 
number of photons is conserved, the isotropic equivalent 
luminosity of the Compton scattered photons is of order 
T 4 times that in the optical/UV photons. If the latter is 
of order the Eddington luminosity for a 10M Q star, then 
the observed isotropic luminosities for bright bursts, 10 13 
to 10 15 L e dd require Lorentz factors exceeding 10 3 . Sha- 
viv and Dar (1995) predicted in their model that inverse 
Compton emission is expected to emerge polarized with a 
correlation between polarization and other quantities asso- 
ciated with the angle of observation. This should include 
overall fluence and peak luminosity, which are down from 
the peak values by factors of 8 and 16 respectively when 
the polarization is at a maximum. The high polarization 
is possible only when the angular width of the beam is 
narrower than 1/L, and this constrains the solid angle of 
the beam to be extremely small if high net polarization is 
to be observed. 

The Compton upscattering model for polarization (Dar 
and Shaviv 1995) should be compared and contrasted to 
the ensheathmcnt model of Eichler and Levinson (1999). 
Both models predicted linear polarization before the an- 
nounced discovery (Coburn and Boggs 2003) In the en- 
sheathment model, the energy of the photon is not neces- 
sarily increased by the scattering and certainly not by a 
dramatic factor. The photons can thus originate in a much 
smaller region without violating the black body limit, and 
short timescales are, in principle, possible even for the 
once scattered radiation (the criticism of Lazzati et al. 
on this point is therefore not justified). The scattering 



material in the enshcathment model may have a relatively 
moderate bulk Lorentz factor and the solid angle of the 
scattered light can subtend a much larger angle than the 
high T scattering material invoked by the Compton up- 
scattering models. It is therefore possible to combine a 
large (lab frame) scattering angle and high polarization. 
It is also easy in the ensheathment model to have a ge- 
ometry in which only scattered radiation is seen by the 
observer. As proposed by Eichler and Levinson (1999), 
the large reduction of the "on-axis" luminosity that was 
apparently achieved by GRB980425/SN1998bw, perhaps 
4 or 5 orders of magnitude, comes with only modest re- 
duction of frequency. It can result from a small fraction of 
photons being scattered, by a non-relativistic or modestly 
relativistic scatterer, into such a large angle as to be ob- 
servable by us. The smooth light curve that was observed 
for GRB980425, rare among GRB's and far removed from 
the relation of Atteia et al. (2003), is consistent with the 
light having reached us via many different paths. 

In recent years, the view that GRB's arc collimatcd (e.g. 
Levinson and Eichler 1993)has been almost universally ac- 
cepted, and the collimating matter is probably a host star 
or a wind. This confirms the assumption that GRB fire- 
balls pass through some sort of surrounding wall in their 
early stages. This has particular significance for models 
in which the early fireball is baryon-free, because it pro- 
vides a source of baryon seeding that could account for 
afterglow-generating plasma escaping intact to infinity. In 
the first place, some baryons are introduced downstream as 
entrained matter or pick-up neutrons that diffuse into the 
outflow from surrounding walls of baryonic matter. This 
matter can be a host star or simply a baryonic wind driven 
from the accretion disk of the black hole. The baryonic 
density in this picture has a complicated transverse struc- 
ture (Levinson and Eichler 2003), the density increases 
outward, and the photosphere can be close to the pair an- 
nihilation radius near the axis, about 10 10 cm, and near 
the walls at the edge (see Fig. 1). 

In this paper, we consider another source of baryons - 
the collimating walls themselves - and note that, due to 
the impact of the GRB photons on the walls, relativistic 
motion of the walls can ensue. This matter can then be 
considered part of the GRB. It can cause highly polarized 
gamma ray emission in certain viewing directions and can 
cause afterglow. It might even be the case that there is a 
smooth transition from the supernova ejecta to the GRB 
fireball. This could account, for example, for the pecu- 
liar features of SN 1998bw, which has a huge expansion 
velocity (c/6) and peculiar abundances. 

This was considered by Eichler (2003, hereafter E03) 
under simplifying geometric assumptions: that the pho- 
tons emerged from a small, effective point source. It was 
argued that the photon pressure would accelerate the wall 
material to a Lorentz factor of 1/ sin\, where \ is the angle 
between the photon flow and the wall, and that the polar- 
ization of scattered photons off the wall would be greater 
than 0.6 with considerable probability. Here we relax the 
point source assumption. In section 2 we consider the 
terminal Lorentz factor in the presence of a photon field 
from a finite source, whose shape is characterized by two 
parameters. In section 3, we consider the effect of finite 
source size on the expected polarization. Finally, we con- 
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sider multiple scattering within a localized neighborhood 
of the point of photon impact on the sheath. 

We demonstrate that in the present model, large polar- 
ization (of order 40 percent) is typical even when the beam 
"acausal", i.e. when it has a much larger opening angle 
than 1/r. This is not the case in most models where the 
intensity along the line of sight typically receives contribu- 
tions from a surrounding neighborhood of velocity vectors 
within the beam, whereby the polarizations of the various 
contributing beams cancel. In the limit of axisymmet- 
ric contribution around the line of sight, this cancellation 
would be complete. In the present model, on the other 
hand, there is a lower degree of symmetry (hence the calcu- 
lations are unfortunately more complicated), because the 
incident radiation arrives from a direction much different 
(even in the comoving frame) from the velocity vector. 

2. POLARIZATION BY A COMPTON SAILING SHEATH 

We consider the dynamics of the inner layers of the 
baryon rich sheath surrounding the baryon poor core, and 
the consequent polarization of the scattered gamma-ray 
emission. We envision that the sheath is irradiated by the 
gamma-ray source which is located at small radii, as il- 
lustrated in Fig 1. The radiative force would accelerate 
the inner layers, leading to a strong shear of the boundary 
layer. Viscous forces may counteract the radiative force, 
and the profile of the Lorentz factor across the boundary 
layer would, in general, depend on the detail of momentum 
transfer between the moving layers of the sheath. In the 
limit of a highly super-Eddington primary photon luminos- 
ity, it seems likely that the radiation pressure overwhelms 
the viscous forces on the exposed layers of scattering mate- 
rial. This case has been considered in E03, and is explored 
further below. 

2.1. Dynamics 
The stress-energy tensor of a magnetized flow is given 

by 

T af3 = hnU a UP - pn afl + -L(F aa F^ + ^^F 2 ), (1) 

4tt 4 

where F^ v is the electromagnetic tensor, U a = (T, T[3) 
is the 4-velocity of the bulk fluid, and n, p, and h, are 
the proper particle density, pressure and specific enthalpy, 
respectively. The energy, momentum and continuity equa- 
tions can be written as 

d a T a » = S", 

d a (nU a ) = 0. (2) 

The above set of equations must be augmented by Maxwell 
equations for the electromagnetic field. The source term 
accounts for all external forces acting on the MHD flow, 
and is expressed here as the sum of a radiative force and all 
other forces (e.g., viscosity mediated by some instabilities, 
electromagnetic stresses, or diffusing neutrons; see Levin- 
son and Eichler 2003): S» = S» + S%. In terms of the dis- 
tribution functions of gamma rays, / 7 , and electrons (we 
don't distinguish here between electrons and positrons), 
/ e , the source term associated with the radiative force is 
given, in the limit of Thomson scattering, by (see e.g., 
Phinney 1982), 



j u j Pj 

(3) 

where p£ and pf^ are the 4 momenta of electrons and 
gamma rays, respectively, as measured in the Lab frame, 
and m m = Pe/ m e is the corresponding electron 4 velocity. 
The gamma-ray distribution function satisfies a Boltzman 
equation: 

P l ^d l _ l f 7 (p 1 ,x) = C 1 {f 1 ,p 1 ,x), (4) 

where the collision operator must satisfy, / pH / C 1 d 3 p 1 = — Sg, 
by virtue of energy and momentum conservation. Once 
5*2 is specified, the above set of equations, subject to ap- 
propriate boundary conditions, can, in principle, be solved 
to yield the structure and dynamics of the boundary layer. 
A complete treatment will be presented elsewhere. 

In what follows, we consider, for illustration, the case 
of a sufficiently cold fluid dominated by radiative forces. 
By sufficiently cold we mean that to a good approxima- 
tion the electron distribution function is given by f e = 
Tn e S 3 (p l e — meU 1 ), where n e is the electron proper den- 
sity. We further suppose that the system is stationary and 
axially symmetric. Projecting eq. (2) on the direction 
perpendicular to the wall, defined by the unit vector n s , 
and using eq. (1), yields the pressure gradient required 
to support the corresponding component of the radiative 
force: 

n s ■ V(p + B 2 /8n) = _ gT ^ r j dv j d£lI v (l-Pn) sin (9. 

(5) 

Here, I v — pi j hcf 1 is the gamma ray intensity, and \i = 
cos 9 denotes the cosine of the angle between the direction 
of impinging photons and the fluid velocity. The zeroth 
component of eqs. (2) yields 

p-V(hT + B 2 /Ann) = -^^x 

nc 

j dv J dQI u {Q){l - fa) [r 2 (l - pn) - 1], (6) 

where the continuity equation has been used. 

The intensity at a given point on the wall generally re- 
ceives contribution from both direct illumination of the 
radiation source and from reflection from other parts of 
the sheath that are exposed to the radiation source. To 
calculate the intensity we must therefore solve a trans- 
fer equation that accounts for multiple scattering between 
the different parts of the wall, as well as for diffusion of 
photons outward into the denser parts of the sheath. Since 
the reflected emission is beamed due to the relativistic mo- 
tion of the sheath, it is clear that eq. (6) and the transfer 
equation obeyed by the intensity are coupled. The solution 
of the radiative transfer equation is beyond the scope of 
this paper. To simplify the analysis, we consider, in what 
follows, direct irradiation of the wall by an effective radi- 
ation source of a finite extent. From eq. (6) it is evident 
that for sufficiently large incidence angles of the impinging 
photons, the fluid will quickly accelerates to its coasting 
speed at which the component of the radiative force in the 
direction of motion nearly vanishes (E03). 

In order to calculate the Lorentz factor of the coasting 
layer, we suppose that the incident radiation is emitted 
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from a surface defined by S — S(x,y,z). The coordinate 
system adopted is shown in Fig. f. The z axis is cho- 
sen to coincide with the symmetry axis. We denote by rp 
the position vector of a fluid element, and by r 7 the posi- 
tion vector of a point on the gamma-ray emitting surface. 
The direction of an incident photon emitted from a point 
defined by r 7 is then given by (see Fig. 1) 



fluid element it is given by /^(r^, k[, v 1 ) = 6fl in (r 7 , h, v) 



(7) 



The geometry of the radiation source adopted in the cal- 
culations presented below is shown in Fig. 2. The gamma- 
ray emitting surface has a geometry of a ring with inner 
and outer radii denoted by a\ and respectively. The 
ring lies in the (x,y) plane at z = 0. The radius of the 
wall at z = is denoted by d, and is normalized such that 
d = 1. The baryonic wall is taken to be conical for conve- 
nience. We note, however, that the difference between a 
conical wall and a curved one (as in the case of collimation) 
is only in the definition of d. In the latter case d represents 
the radius of a fictitious cone that extends from the loca- 
tion of the scattering material along the velocity vectors. 
The intensity in the source is taken to be homogeneous, 
that is, independent of x and y inside the ring (and zero 
outside, of course). The spectrum of incident radiation is 
taken to be a power law with index a — 0.5 (Jj„ oc v~ a ). 
The ring emission is allowed to be beamed. The beaming 
is parametrized by an angle 6 1 , measured with respect to 
the symmetry (z) axis. For a given beaming angle 7 , only 

rays that satisfy the condition ki • z — cos(# 7 ) for the ki 
given by eq. (7) can scatter off a fluid element located at 

In the coasting regime the R.H.S of eq. (6) nearly van- 
ishes, and the solution of eq. (6) reduces to 

= n-( K 2 -l)^ 2 , 



with 



K = 



J dv J dn s i v (i + fj, 2 ) 

2 J dv J dO-s^Iv 



(8) 



(9) 



where the inner integration is carried over the gamma-ray 
emitting surface, S, with dfl s — (dS ■ rg/r^). 

In the limit a\ = a 2 = that corresponds to a point 
source located on the symmetry axis, the intensity is given 
by I oc S(p — fi ), with /i = z/(z 2 + rf 2 ) 1 / 2 , and eq. (8) 
reduces to = no, recovering the result obtained by E03. 

2.2. Polarization 

Consider now the polarization of the radiation scat- 
tered by the moving wall. Using the coordinate system 
depicted in Fig.l, we express the velocity vector of a 
fluid clement, measured with respect to the Lab frame, as 
= 0(sin 9g cos(j)fjX+sm8p sin <fi gy + cos 6 pz), and the line 
of sight direction as h = sm9 n cos<t> n x + sin #„ sin + 
cos 9 n z. The vector fcj given by eq. (7), and the vector n 
are transformed into the comoving frame of the fluid cle- 
ment (henceforth denoted by prime), using the aberration 
formula: cos (9^ = k[ ■ = {h ■ - 0)/[l - 0{h ■ 0)], and 
likewise for cos 9' n3 = n' ■ 0. 

Let Iin(ry, h, v), with ki given by eq. (7), denote the 
intensity of the incident radiation. In the rest frame of the 



where Si = v'/v = [T (1 + cos 0' i/3 )] 1 is the corresponding 
Doppler factor. The polarization vector in the comoving 
frame lies along the direction t = h' x ki /\n' x fcj |. In 
terms of the unit vector e = n' x 0/\n' x 0\, the Stokes 
angle is cosx = t ■ e. The comoving emissivity is given by 



j'p = -^ a T< J iL(^ + cos^' 2 )dn s , 

and the corresponding Stokes parameters by 



(10) 



9 

u 



—a T n' e [ dn s l' m (l~cosiP' 2 ) ( C0S l X 
16tt e J my ' \ sm2x 



(11) 

where cos ip' — n' ■ ki . If we now measure the Stokes angle 
with respect to the fixed vector b = n x z/\h x z\, then 
the Stokes parameters are rotated by an angle — 2rj, where 
cos Tj = b ■ e, viz., 
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q'p cos 2rj — u'p sin 2r/ 
q'p sin 2rj + u'n cos 2rj 



(12) 



Transforming to the Lab frame and averaging over direc- 
tions of emitting fluid elements, we finally obtain 



I = 



16tt 



fff J^nd^dOp, (13) 



and 



U 



3 
16tt 



q' g cos 2r) - u' g sin 2r] \ k 

(14) 

where r = arn e l is the optical depth along the sight line, 
and 8 n = [T(l — cosOnp)]" 1 is the corresponding Doppler 
factor. The factor 1/r is introduced because of the trans- 
formation of the electron density from the rest frame of the 
fluid element to the Lab frame. Here we allow for a spread 
in the opening angle of the conical wall. The dependence 
of n e and T on the angle 9g should be specified in gen- 
eral. Below we examine two cases: a thin cone, defined by 
n e (0p) oc S(9fj — 9q), and a cone with a constant density in 
the region 9 — A9 < 9p < 9o + A9 and zero outside this 
range. The index k depends, in general, on the spectrum 
of incident radiation, as well as on the kinematics of the 
emission region, and is expected to lie in the range from 
2.5 to 3.5. The dependence of the polarization on k is ex- 
plored below. Finally, the polarization degree is given by 
P = (Q 2 + U 2 ) 1 / 2 /!. We note that for an axially sym- 
metric distribution of the scattering material, the Stokes 
parameter U should vanish (Bcgclman & Sikora, 1987). 
As a check, we numerically computed U in each run and 
verified that, within the accuracy of our computation, it 
is zero. 

We also calculate below the probability for observing a 
source with a polarization P > tt, given explicitly by, 



x(P > tt) 



2tt 
~A 



r3/2 



dfx, 



P(li)>n 



where the integration is carried over all solid angles for 
which the polarization fraction P exceeds the value tt, and 
A = f P r \ >0 Vmaxd^l is a normalization constant. 
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3. RESULTS 

Eqs. (10) through (14) have been integrated numeri- 
cally to yield the intensity and polarization of the radia- 
tion scattered off a section of the wall that moves with a 
given velocity. In cases corresponding to Compton sailing 
material, the velocity of the scattering matter was first 
calculated using eqs. (9) and (8). As a check on our code, 
we computed the scattered intensity for a point radiation 
source located on the symmetry axis ( a\ — a 2 = 0), and 
compared the results with the analytic approximation de- 
rived by E03. An example is shown in Fig. 3 for z = 100, 
k = 3, and 9q = 0.1. The corresponding Lorentz factor of 
the Compton sailing fluid is T = 101. The numerical result 
is given by the dashed line and the analytic one by the solid 
line. We also plotted the polarization fraction P (dotted 
line). As seen, the agreement is good. The slight asym- 
metry of the dashed curve around the direction of fluid 
motion, A = 0, is due to second order effects associated 
with the curvature of the wall which are not accounted for 
by the analysis of E03. Note also that the viewing angle 
along which the intensity peaks is slightly shifted from the 
direction A = 0, and in particular does not coincide with 
the direction of maximum polarization. This is due to the 
angular dependence of the Thomson cross section. This is 
the reason for the relatively rapid decline of the x{P > it) 
curve. As a second test case, we calculated the polariza- 
tion for a radiating ring of radii a\ — > a 2 = 1, with a highly 
beamed emission (# 7 << 1). This case corresponds to the 
head on approximation considered in Eichler and Levin- 
son (2003). The Lorentz factor of the moving wall in this 
case is arbitrarily chosen (not the Compton sailing value, 
of course) . The resultant polarization curves are displayed 
in Fig. 4. We find perfect agreement as # 7 — ► 0. 

The dependence of the polarization produced by a 
Compton sailing sheath on the various parameters is ex- 
amined in Figs. 5-7. The upper panel in each figure 
exhibits the polarization fraction P versus the line of sight 
angle 9 n , and the lower panel the corresponding Euclidean 
probabilities x(P > it), given by eq. (15), against tt. 

Note that in the limit that sources are so bright that 
they can be seen anywhere in the universe at most view- 
ing angles, then the probability of a given polarization is 
indicated by A9, i.e. the x axis in the upper panels of 
figures 5-7, except for effects from cylindrical geometry, 
which are small when 9 T » 1 and 9 O /A0 3> 1. In the op- 
posite limit, probably not the most relevant for long, hard 
bursts, we assume Euclidean geometry with V/V max = .5. 
The probabilities for observing polarization 7r almost cer- 
tainly lies between these two limits. 

Fig. 5 shows the dependence of the polarization curves 
on the opening angle of the flow, 9q — (3 ■ z. In this exam- 
ple a 2 = d = 1, that is, the radiating ring extends from 
the wall, and oi = 0.9, as indicated. The ring is assumed 
to radiate isotropically in this example (more precisely, 
9 1 > 2#o is assumed). As seen, the polarization approaches 
zero on the axis (9 n = 0), as required by the axial sym- 
metry. The Lorentz factor of the Compton sailing fluid 
obtained in this run is T ~ 56 for 9 = and , T ~ 57 
for 6>o = 0.1. The reason that T increases slightly with 
increasing 9 n , is that the angle between the directions of 
incident photons and slightly decreases on the average, 
but since the ring radiates isotropically and the distance z 



from the ring to the scattering material is rather large, the 
effect is small. The polarization peaks at 9 n = 9q (but is 
not 100% owing to the contribution from elements moving 
at different (j>p), as long as 9q is not too small. Quite gener- 
ally we find that the maximum polarization is suppressed 
for sufficiently small 9 , as seen in Fig. 5, by virtue of 
the axial symmetry. The kinks in the polarizations curves 
correspond to a 90° change in the polarization vector. 

The dependence of the polarization on the index k is 
displayed in Fig. 6. The Lorentz factor obtained for all 
cases shown in this example is 57. As seen, the P curves 
depend weakly on k, and mainly in the side lobes, owing 
to the large beaming factor of the scattered emission. The 
primary reason is that the main bump of the polarization 
curves in all cases lies in the range TA << 1. For much 
smaller Lorentz factors we find somewhat stronger depen- 
dence on k. The intensity is more sensitive to fc, and this 
is reflected in the V max curves, as seen. 

Fig. 7 exhibits the dependence of the polarization on 
the geometry of the incident radiation source. Three cases 
are examined: a point source (a\ = a 2 = 0), a disk extend- 
ing from the wall (ai = 0, a 2 — 1), and an infinitely thin 
ring (ai = a 2 = 1). As seen, a point source gives rise to 
a nearly maximum polarization in the direction 9 n = 9q, 
as expected in the limit of large beaming. This is due to 
the fact that all incident rays are perpendicular to the di- 
rection of motion in the comoving frame. The peak of the 
polarization curve is smaller in the case of an extended 
source by virtue of the finite range of incidence angles, 
but is still rather high for sufficiently large Lorentz factors 
(r = 64 in the case of the disk and 56 for the thin ring). 

In Fig. 8 we explore how the polarization depends on 
the Lorentz factor of the Compton sailing matter. All 
three cases shown were calculated using a point radiation 
source (a\ = a 2 = 0). The source is located at a dif- 
ferent distance from the scattering material in each case. 
The Lorentz factors of the sailing material, roughly equal 
r = (z 2 + d 2 ) 1 / 2 /d are indicated. As expected, the maxi- 
mum polarization is reduced for small Lorentz factors, ow- 
ing to the larger contribution from fluid elements located 
at larger angular separation with respect to the sight line. 
We found, however, that the cumulative probability curves 
converge to the solid line for T > 25, as can be inferred 
from the figure. 

Fig. 9 examines the effect of Compton sailing. Here 
we compare the cumulative V max curve obtained in the 
Compton sailing case (indicated as T c in the figure) with 
cases obtained when the Lorentz factor deviates from this 
value. We also plotted for comparison the result obtained 
in the head on approximation in the limit T9o — ► (Eich- 
ler and Levinson 2003). It is seen that Compton sailing 
greatly increases the likelihood of strong polarization. 

The results described above correspond to an infinitely 
thin cone with an opening angle 9q. The effect of a spread 
in the opening angle of the cone of optically thin scattering 
material is examined in Fig. 10. In this example the veloc- 
ity vectors are assumed to lie in the range between 9q — A9 
and 6*o + A9. All fluid elements are assumed to move with 
their Compton sailing Lorentz factor. As seen, the polar- 
ization is lowered when A9 is much larger than 1/r, but, 
in the absence of any cylindrical symmetry around the line 
of sight, does not plummet to zero in this limit. 
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Finally, the effect of photon diffusion in cases of a large 
optical depth is depicted in Fig. 11. Solving the radia- 
tive transfer in a medium with a modest optical depth is 
a formidable task, and is beyond the scope of this paper. 
Here we consider diffuse reflection from a plane parallel, 
semi-infinite slab. We suppose that in the comoving frame 
the radiation incident normal to the wall, and use the re- 
sults derived in Chandrashekar (1960) for the polarized 
intensity reflected off the wall. We then transform the in- 
tensities back to the Lab frame in order to compute the 
probability x(P > it). The upper panel exhibits the polar- 
ization curve and the (normalized) Euclidean V max curve, 
and the lower panel the cumulative probability as usual. 
The direction of maximum polarization is along the wall 
of course, with peak polarization of 0.47. The total in- 
tensity of the reflected radiation increases with increasing 
angle from the wall, and is maximum at an angle of 1/T, 
as measured in the Lab frame (about a factor of 2 larger 
than along the wall) . This is reflected in the shape of the 
V m ax curve. 

4. CONCLUSIONS 

We have shown, under less simplified assumptions than 
made in Eichlcr (2003), that when matter is Compton 
dragged into a "sailing" equilibrium by photons, the po- 
larization can be quite high for singly scattered photons. 
The direction of maximum polarization coincides with the 
maximum intensity of the polarized component. Finite 
source size is found to not qualitatively alter this conclu- 
sion. Finite source angular width 58, which can nearly 
eliminate net polarization in most models, is found in the 
present model to reduce it only somewhat, and to leave 
about 40 percent polarization for most relevant viewing 
angles. 

In the picture presented here, singly scattered photons 
may in principle compete with unscattered ones, because 
the photon energy as seen in the lab frame is not dra- 
matically altered. [This is in contrast to the scenario of 
Shaviv and Dar, 1995, Lazzati et al. 2000, where the pho- 
ton energy is dramatically raised by the scattering, so that 
virtually all the gamma rays are attributed to such scat- 
tering.] However, the scattered and unscattered photons 
do not necessarily arrive at the same time, are not con- 
centrated in the same directions, and do not necessarily 



have exactly the same spectrum. [In some cases, the line 
of sight does not necessarily admit unscattered photons, 
e.g. as was argued (Nakamura 1998, Eichler and Levinson 
1999) for the case of GRB980425, however, this GRB may 
be exceptional.] It is thus important that future GRB po- 
larization experiments be sensitive to the polarization as 
a function of time and energy range. 

If the photons are scattered (undergoing " diffuse reflec- 
tion" in the terminology of Chandrasekhar, 1960) off an 
optically thick sheath, then the maximum polarization cor- 
responds to 47 percent along the direction of maximum 
intensity of the polarized component. Typically, polariza- 
tions of 30 to 47 percent are expected in this situation 
among the scattered photons. Consistently observed val- 
ues in this range of polarizations, with a sharp cutoff at 
47 percent, might thus be a signature of diffuse reflection. 

The situation in which the optical depth of the scatter- 
ers is slightly greater than unity remains to be studied in 
detail. We believe that the results will be not too much dif- 
ferent from the case of single scattering, because the path 
of maximal polarization corresponds to minimum optical 
depth in a relativistically expanding flow. In any case, 
they should lie between the two extremes considered in 
this paper. It might be claimed that an optical depth of 
order unity is "fine tuned". However, in the present pic- 
ture, the sheath tapers to infinity, so most of its area is in 
fact comprised of material with an optical depth of order 
unity. The path of least optical depth is along the direc- 
tion of motion, because then the relative velocity between 
the photon and the material, as seen by the lab frame is 
the smallest (Levinson and Eichlcr, 2003). 

If future polarization studies eventually show some indi- 
cation of scattering of GRB photons, correlation between 
maximum polarization and maximum intensity might pro- 
vide clues as to whether energy is indeed flowing from 
the scattering material to the photons or whether most of 
the scattering occurs after Compton equilibrium has been 
reached, which would suggest that matter is being dragged 
by the photons rather than upscattering them. 
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tional Science Foudnation, the Arnow Chair of Theoreti- 
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grant for a Israeli Center for High Energy Astrophysics. 
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Fig. 1. — Schematic illustration of scattering by a relativistic sheath. Gamma rays are emitted from a surface marked as S in the figure, and 
scatter off matter that ensheath the ultrarelativistic, optically thin core. The symmetry axis of the system is taken to be in the z-direction. 
The position of a scattering center on the wall, moving with a velocity /3, is defined by the vector rp. The position vector of a point on the 
gamma-ray emitting surface is labeled by r 7 . The sight line direction is denoted by h. 
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Fig. 2. — A sketch of the radiation source geometry adopted for the calculations is shown. The emitting surface has a geometry of a ring 
with inner and outer radii denoted by ai and 02, respectively. The ring lies in the (x, y) plane at z = and is centered around the symmetry 
axis of the system. The radius of the wall containing the scattering material at 2 = is denoted by d. 




Fig. 3. — Scattered intensity (dashed lines) and polarization degree (dotted line), calculated using a point radiation source located on the 
symmetry axis, are plotted as a function of viewing angle 8 n . A plot of the intensity calculated using the analytic approximation derived in 
Eichler (2003) is shown for a comparison (solid line). The parameters 8q and k are defined in the text. 




Fig. 4. — The polarization degree as a function of T(sin 9q — sin 9), calculated numerically using a radiating ring with a\ = C12 = 1 and highly 
beamed emission for different values of T9q, is shown for the purpose of comparison with the results obtained in the head on approximation 
by Eichler and Levinson (2003). 
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Fig. 5. — The upper panel exhibits polarization curves for different values of the opening angle 6q. The polarization in the inner bump is 
parallel to the wall, and changes by 90° across the kinks. The bottom panel displays the corresponding probability for observing a source 
with polarization higher than w, as given by eq. (15) in the text. 
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Fig. 7. — Same as Fig. 5, but for different ring parameters 




Fig. 9. — Comparison of the cumulative probability curves obtained for a Compton sailing wall (indicated by r c ), with those obtained in 
cases where the Lorentz factor of the scattering material deviates from the Compton sailing value (see text for further discussion). 




Fig. 10. — The polarization produced by scattering off a matter moving on conical shells encompassing a range of opening angles between 
00 — A9 and 9q + A8. Each sub-shell is coasting with a Lorentz factor corresponding to the Compton sailing value. The Loerntz factor of the 
central sub-shell is indicated.. 
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Fig. 11. — Effect of multiple scattering inside the wall on the polarization of emergent radiation. Shown are the polarization of a radiation 
reflected off a semi-infinite, plane-parallel slab, as a function of the angle A from the wall (upper panel), and the cumulative probability (lower 
panel) . The direction of the incident beam is taken to be normal to wall in the wall's rest frame. 



